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The superconductor to normal metal phase transition induced by the in-plane magnetic field is 
studied in free-standing Pb(lll) nanofilms. In the considered structures the energy quantization 
induced by the confinement leads to the thickness-dependent oscillations of the critical field (the so- 
called ’shape resonances’). In this paper we examine the influence of the orbital effect on the in-plane 
critical magnetic field in nanofilms. We demonstrate that the orbital term suppresses the critical 
field and reduces the amplitude of the thickness-dependent critical field oscillations. Moreover, due 
to the orbital effect, the slope ~ Tc at Tc(0) becomes finite and decreases with increasing film 
thickness in agreement with recent experiments. The temperature t* at which the superconductor 
to normal metal phase transition becomes of the first order is also analyzed. 


I. INTRODUCTION 


The huge progress in nanotechnology which has been 
made in the last decade reopens the issue of super¬ 
conducting properties of metallic nanostructures i.e. 
nanofilmsi^— , nanowires^^— or metallic grains^rii. The 
studies of the quantum size effect and its influence on 
the paired phase in thin films was initiated by Blatt and 
Thomson in 1963^^. The main finding of their worbi^ 
was the appearance of sharp oscillations of the critical 
temperature as a function of the film thickness. As it 
was arguedi^ this effect results from the confinement of 
the electron motion in the direction perpendicular to the 
film. If the size of the system becomes comparable to the 
electron wave length, the Fermi sphere splits into a set 
of discrete two-dimensional subbands, energy of which 
increases with decreasing film thickness. Each time a 
bottom of a subband passes through the Fermi level, 
a sharp peak of the critical temperature appears. Due 
to technological difficulties in the preparation of uniform 
films, which were typically polycrystalline and contained 
a large number of defects, the experimental observation 
of the so-called shape resonances has been reported only 
recentlyi^^— . Measurements of the critical temperature 
for Pb(lll) nanofilms grown on Si substrat o^^d^ revealed 
that the period of the thickness-dependent oscillations is 
equal to ^ 2 ML. Additionally, the beating effect with the 
periodicity varying from 7 ML to 14 ML was observed. 
This feature, called bilayer or even-odd oscillations, was 
theoretically studied by Shanenko et al. in Ref/^. The 
studies of superconducting properties of the ultra-thin 
films have been recently extended to the case of high 
temperature and multiband superconductors. Recently, 
the enhancement of the superconducting critical temper¬ 
ature with respect to the bulk limit has been reported for 
interfaces and heterostuctures based on cupratesi^, iron 
pnictides— and LaAlOa/SrTiOg^. Moreover, the exper¬ 
imental reports on the growth techniques of high quality 


MgB 2 films of thicknesses less than 10 nm—>22 entailed 
a series of theoretical papers describing the quantum size 
effect in multiband thin film superconductors—^—. 

Recent studies devoted to superconductivity in the 
nanoscale regime concern the effect of the quantum 
confinement on the superconductor to normal metal 
phase transition induced by the magnetic field^^. The 
thickness-dependent oscillations of the perpendicular and 
parallel critical field for ultra-thin lead films were re¬ 
ported by Bao et al. in Ref—. Moreover, the study of the 
superconductor to normal metal transition induced by 
the parallel magnetic field for Pb monolayer was recently 
presented by Sekihara et al. in Ref—. In both experi¬ 
ments, the measured parallel critical field was higher than 
the Pauli paramagnetic limit. This unusual behavior has 
been explained in our recent paper s^^i^^ in which we have 
investigated the quantum size effect on the in-plane crit¬ 
ical field in paramagnetic limit. We have shown^ that 
the zero-temperature critical field for nanofilms is higher 
than the Clogston - Chandrasekhar (CC) paramagnetic 
limit and diverges to the CC (Pauli) limit for sufficiently 
thick films. This fact has been explained on the basis of 
the spatially varying energy gap induced by the confine¬ 
ment. In Ref— the new formula for the paramagnetic 
critical field in nanofilms has been proposed. However, 
the analysis presented in Ref— has been carried out in 
the paramagnetic limit. The extension of this study for 
films thicker than 15 ML requires the inclusion of the or¬ 
bital effect which significantly affects the superconductor 
to normal metal transition. According to our knowledge 
such study has not been reported until now. 

In the present paper we consider free-standing Pb(lll) 
metallic nanofilms and investigate the orbital effect on 
the superconductor-normal metal transition driven by 
the in-plane magnetic field. Based on the analysis of the 
spatially dependent energy gap we study the influence of 
the orbital effect on the critical magnetic field oscillations 
induced by the confinement. The analysis of the thermal 
effect in terms of the orbital effect is also included. The 
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paper is organized as follows: in Sec. |TT] we introduce 
the basic concepts of the theoretical scheme based on the 
BCS theory, in Sec. IIIII we present the results while the 
summary is included in Sec. lYl 


II. THEORETICAL METHOD 


The phonon-mediated superconductivity in metallic 
nanofilms can be described with the use of the BCS the¬ 
ory. The Hamiltonian of the system is given by 


a 
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( 1 ) 


where a indexes the spin state (f, i) and g is the electron- 
phonon coupling. In the presence of the in-plane mag¬ 
netic field id||, the single-electron Hamiltonian can 
be expressed as 


+ “A.) -b s/rsidii — (2) 

where s = -bl(—1) for a =t (i), m is the effective elec¬ 
tron mass, fj,F is the chemical potential, A = (0, —id||Z, 0) 
is the vector potential corresponding to the magnetic field 
H = (7d||,0 , 0) applied in-plane and the energy gap A(r) 
is defined as 


A(r) =- 5 <^^(r, 4 .)^i(r,t)) . ( 3 ) 

In ultrathin nanofilms the electron motion is limited in 
the direction perpendicular to the film (z axis) resulting 
in the quantization of the electron energy. We assume 
that the system is infinite in the x — y plane. Thus, the 
field operators in Eq. o are expressed as 

^ (r, cr) = ^ (/)kn (r) Ckna, ( 4 ) 

n,k 

cL<t^ (5) 

k.n 

where Ckrtcr(Ck„c,) anihilation (creation) operator 

for an electron with spin u in a state characterized by 
the quantum numbers (k, n) while ^kn (r) is the single¬ 
electron eigenfunction of the Hamiltonian idf whose ex¬ 
plicit form is given by 

0kn(r) = (6) 

where k = {kx,ky) is the electron wave vector and n 
labels the discrete quantum states induced by the con¬ 
finement along the z axis. 

By using the eigenfunctions given by Eq. (jS]), one can 


reduce the Hamiltonian ([^l to the ID form which corre¬ 
sponds to the z dependent part - 


"4^ 2m dz^ 
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Note, that in the presence of the magnetic field id|| the 
eigenfunctions (fk^niz) depend on the ky component of 
the wave vector. In our calculations ipk^niz) are deter¬ 
mined numerically by the diagonalization of the Hamil¬ 
tonian © in the basis of the quantum well states 


‘Pkyniz) 
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where d is the film thickness whereby we adopt the hard- 
wall confinement as the boundary condition in the z di¬ 
rection. 

By using the Bogoliubov - de Gennes transformation 
Ckn<7 = Ukna'lkn + the energy gap in the band 

n defined as A„ = (())kn|A(r)|(?ikn) can be expressed as 
follows 

An' = j dkxdky X (9) 

E gfe.n>n .-f" [1 - nKu) - nKu)] > 

n' Skn + 

where ^kn is the single-electron energy, = 

A\/?kra + |Anp is the quasi-particle energy, f{E) is the 
Fermi-Dirac distribution and Ckyu'n are the interaction- 
matrix elements given by 


Ckyn'n = / dz(pkyn'{z)(p-kyn'{z)(pkyn{z)(p-kyn{z). 

( 10 ) 

The summation in Eq. (US is carried out only over the 
single-electron states with energy ^kn inside the Debye 
window l^kral < where wf is the Debye frequency^. 
The spatially dependent order parameter A(z) can be 
expressed as 


A('2) = ^ / dkxdky'^ipkyniz)(f-kyn{z) X 

[^-fiKn)-fiKn)]- ( 11 ) 


A, 


2\/?L + A^ 


The energy gap A(z) for given nanofilm thickness d is 
calculated in the self-consistent manner using the fol¬ 
lowing procedure. For each value of the wave vector ky 
(note that the range of the wave vector is limited by the 
condition |<Ckn,| < huJo) we calculate the single-electron 
wave functions ipkyn{z) and the energies ^kn by the diag¬ 
onalization of the Hamiltonian ([2]) in the basis given by 
Eq. ([5]). Then, the single-electron wave functions (/Sj, „(z) 
are used to determine the staring value of 5n- The final 
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value of 6n is calculated by the self-consistent procedure 
given by Eq. [TUI Finally, having (5„ and the spa¬ 

tially dependent order parameter A(z) is calculated by 
the use of Eq. |TT] The calculations in the paramagnetic 
limit (with no orbital effect) are carried out by putting 
iT|| = 0 in Eq. |7| everywhere except the Zeeman term. 
In this case the sum (|5]) reduces to a single term - the 
situation corresponds to a infinite quantum well with the 
spin Zeeman splitting. 

We should note that aforementioned procedure can 
lead to solutions with A 7 ^ 0 even for the values of 
the magnetic field for which the free energy of the 
superconducting state is greater than the free energy 
corresponding to the normal metal solution (A = 0). 
Therefore, for each value of the magnetic field we 
calculate and compare the free energies of the normal 
and superconducting phase. The theoretical study of 
the free energy in superconducting nanostructures has 
been presented in detail in Refi^. 

Due to the confinement the chemical potential can 
strongly deviate from the bulk value. For this reason 
for each nanofilm thickness we determine the chemical 
potential from the formula 

— ~d j' ^ ^ j' l^knfT^fcyn('2^)| f 

+ \Vkna‘fkyniz)f[l - f{Ei,n)]^, ( 12 ) 

where Ue is the electron density corresponding to the bulk 
value (corresponding to the chemical potential Hbuik)- 


III. RESULTS 

In the present paper we consider the free-standing 
Pb(lll) nanofilms. The first-principle calculations 
of the quantized band structure for Pb nanofilms in 
( 111 ) and ( 100 ) directions are presented in Refsj^^— . 
Authors of these papers have pointed out that in (111) 
direction the energy dispersion is nearly parabolic 
and the quantum size effect can be well described by 
the quantum well states centered at the L-point of 
a two-dimensional Brillouin zone^. Based on these 
results, in our analysis we use the parabolic band 
approximation treating the bulk Fermi level Hbuik and 
the electron mass m as the fitting parameters. Their 
values are determined based on the results from the 
first-principle calculations for Pb(lll) presented in 
Refs.— 1 ^. We use the following values of the parame¬ 
ters: gNbuik{0) = 0.39 where NbuikiO) = mkr/i^Tr'^h'^) 
is the bulk density of the single-electron states at the 
Fermi level, Tlcod = 8.27 meV, the bulk energy gap 
Ab„;fc = 1.3 meV and g,buik = 3.8 eV which corresponds 
to the electron density rif. = 4.2 x 10^^ cm“^. 


A. Orbital effect on the critical magnetic field. 

Figure [T] displays the in-plane critical field i7c,|| in 
units of Hpf as a function of the nanofilm thickness d, 
where is the paramagnetic critical field calculated 
on the basis of the Clogston-Chandrasekhar limit^i2S 
= Abuik/V^gB) which for bulk Pb gives = 

15.9 T. The thickness range under consideration is chosen 
on the basis of the experiments which report the stable 
Pb(lll) nanofilms with the thickness varying from 5 ML 
to 30 Within our analysis the value of i7c,|| for 

each nanofilm thickness is defined as the field for which 
the spatially averaged energy gap A = (l/d) A(z)dz 
drops below O.OlAbuife. The physical origin of ’the tooth- 




FIG. 1. (Color online) In-plane critical magnetic field Rt:,|| 
in units of as a function of the nanofilm thickness d 

calculated with and without (Pauli approximation) the in¬ 
clusion of the orbital effect. Results for the temperature 
(a) T = 0 K and (b) T = 6 K. Insets present i7t,,|i(d) as 
a functions of number of Pb monolayers assuming the lattice 
constant apb = 0.286 nm. Green squares in panel (a) denote 
the film thicknesses chosen for presentation in Fig. (5) 

like’ oscillations in Fig.|T]can be explained in terms of the 
electron energy quantization related to the confinement 
of the electron motion in the direction perpendicular to 
the film^i2^. If the thickness of the nanofilm becomes 
comparable with the electrons wave length the Fermi 
sphare splits into the set of discrete subbands. Increase 
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of the film thickness results in a decrease of the discrete 
subband energies. Each time a subband passes through 
the Fermi level the density of states in the energy win¬ 
dow [^ — fiojD, d]^ in which the phonon-mediated 

pairing occurs, abruptly increases. This phenomenon 
leads to the thickness-dependent oscillations of the en¬ 
ergy gap and the critical field presented in Fig. [T] As we 
can see, in the thickness range between two subsequent 
resonances, the critical magnetic field almost linearly de¬ 
creases with increasing nanofilm thicknes a^^’^^ . The pre¬ 
dicted enhancement of the critical field reaches almost 
twice the value of the paramagnetic limit in the bulk Pb. 
Similar behavior was recently reported for Pb nanofllms 
in Ref.— where it is shown that experimentally measured 
critical field Hc^\\ was much higher than the paramag¬ 
netic limit Hp^. The increase of the paramagnetic in¬ 
plane critical field in nanofllms has been explained in our 
recent paper—. In the insets of Fig. [I] we show the value 
of ^^c,|| Eis a function of monolayer number (we assume 
the lattice constant apb = 0.286 nm corresponding to 
the bulk value). Results presented in such manner reveal 
the bilayer (even-odd) oscillations with the beating ef¬ 
fect observed in experiments with Pb nanofilm s^^d^ and 
explained in Refs— 

In order to analyze the influence of the orbital effect on 
the in-plane critical field, in Fig. [T] we present the results 
of the calculations carried out with the inclusion of the 
orbital effect (solid, red line) and in the Pauli approxi¬ 
mation (dashed, blue line). From Fig.[T]one can see that 
the orbital effect leads to a decrease of the critical field 
what is clearly visible for thick nanofilms. In compari¬ 
son with the Pauli approximation the amplitude of the 
critical field oscillations is also reduced with increasing 
nanofilm thickness. One should note (see Figl2]) that the 
orbital effect significantly affects the value of i7c,||(d) for 
nanofilms with the thickness greater than 15 ML. The up- 



d (ML) 


FIG. 2. (Color online) Difference between the critical mag¬ 
netic field calculated in the Pauli approximation (Ri(||) and 
with the inclusion of the orbital effect (Rc,||) in units of Hp^. 
Results for T = 0 K. 


per limit of the thickness above which the Pauli approxi¬ 
mation is no longer satisfied can be approximated by the 
magnetic length an = where is the para¬ 

magnetic (Pauli) critical field. Nevertheless, the use of 
the Clogston-Chandrasekhar paramagnetic field for the 
bulk Hp^ = 15.9 T gives an ~ 22 ML which is greater 
that the thickness limit 15 ML determined from the nu¬ 
merical calculations. This discrepancy results from the 
different values of the paramagnetic critical field in 
nanofllms in comparison with the bulk valued. 

The suppression of the critical field induced by the or¬ 
bital effect can be explained on the basis of the classical 
Lorentz force acting on electrons in the magnetic field. 
The parallel magnetic field applied along the x axis (for 




FIG. 3. (Golor online) (a) Wave functions for the first 
and the second ip 2 ,±k[z) quantum well states cal¬ 
culated for opposite k vectors. States with opposite k are 
shifted towards the opposite edges of the sample. The wave 
function ipi^pk{z) = ip\^-k{z) obtained for the case with no 
orbital effect is marked by the black, dotted line, (b) Spa¬ 
tial dependent energy gap A{z) in units of Abuik calculated 
for the magnetic field H\^/Hp^ = 0.7 with (solid line) and 
without (dashed line) orbital effect. Results for d — 30 ML. 


the s-wave superconductors the direction of iL|| is not rel¬ 
evant) results in the Lorentz force directed perpendicular 
to the plane (along the 2 axis). For electrons with oppo¬ 
site k vector (which form the Cooper pair) the Lorentz 
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forces have opposite orientations. As a consequence, for 
nanofilms with the thickness reduced to several monolay¬ 
ers, the electron density in the z-direction is shifted from 
the center to the edges (in accordance with the Lorentz 
force orientation). The shift of the density for electrons 
with opposite k towards opposite edges of the sample 
is clearly visible in Fig. [31(a) which presents the single¬ 
electron wave functions for the first 'fi,±k{z) and the sec¬ 
ond ^P 2 ,±k[z) quantum well state. This effect does not 
take place for the case with no orbital term included for 
which (/?i__i_fc( 2 ) = ipi^-k{z)- According to Eq. m the en¬ 
ergy gap depends on the overlap between the wave func¬ 
tions of electrons with opposite momenta, which from 
the Cooper pairs. As presented in Fig. jSKa) the orbital 
effect results in the reduction of this overlap and conse¬ 
quently leads to the suppression of the energy gap and 
the critical magnetic field. We should mention that the 
presented explanation is correct only for nanofilms which 
are too thin for the formation of the vortex state, which 
is the case considered here. The Lorentz force manifest 
itself also in the spatial distribution of the energy gap 
presented in Fig.jHKb). It is well known that in supercon¬ 
ducting nanostructures the energy gap is not uniform, as 
in the bulk, but it depends on the position. The spa¬ 
tially dependent energy gap for the nanofilm thickness 
d = 30 ML and the magnetic field H\\/Hp^ = 0.7 is 
presented in Fig. |3Kb). The comparison of the results 
calculated with and without the inclusion of the orbital 
effect allows to conclude that the reduction of the gap 
parameter caused by the orbital effect is most significant 
at the center of the film. 



FIG. 4. Energy gap A„ in units of Abuik as a function of the 
magnetic field 77|| in units of Hp^ for the nanofilm thickness 
d = 30 ML. 


Let us now explain the ’step-like’ dependence of the 
difference iL^|| — 77^,1 1 as a function of the film thickness 
presented in Fig. ^ The thicknesses d for which the 
steps appear correspond to peaks of the critical field seen 
in Fig. [TJ This indicates that the presented behavior 


is directly related to the discrete energy spectrum 
induced by the confinement. Therefore a subband 
passing through the Fermi level (with increasing d) 
not only contributes to the increase of the density of 
state at the Fermi level but also to the orbital effect. 
This orbital contribution from subsequent subbands 
manifests as the steps in Fig. |2j Since the orbital effect 
corresponding to each subband is different, the heights 
of the steps also differ. The different importance of 
the orbital term in different subbands is clearly visible 
in Fig. |4| in which we present the subband energy 
gap A„ as a function of the magnetic field iL||. For 
comparison A„(i7||) in the Pauli limit is also shown. 
The slow decrease of the energy gap with increasing 
magnetic field seen in Fig. |4| results from the orbital 
effect. Note that in the Pauli limit A„ is the same for 
all subbands and does not depend on the magnetic field 
in the superconducting state. This clearly show that the 
orbital effect coming from different subbands n is dif¬ 
ferent resulting in the steps with unequal height in Fig. [5] 


1. Orbital effect on i7t,,||(7c) 

In this subsection we discuss the influence of the or¬ 
bital effect on the superconductor to normal metal phase 
transition induced by the magnetic field in non-zero tem¬ 
perature. In Fig. Owe present the spatially averaged en¬ 
ergy gap as a function of magnetic field and temperature 
for the nanofilm thickness (a) d = 21 ML and (b) 23 ML. 
The chosen thicknesses correspond to the maximum and 
the minimum of the zero-temperature critical field pre¬ 
sented in Fig.ITJa). As one can see the range of magnetic 
field and temperature in which the film remains in the su¬ 
perconducting state varies with its thickness. Since the 
significant impact of the orbital effect can be observed 
in the vicinity of T = Tc(0), we restrict our analysis 
to this range. Fig. 0 presents the h — t phase diagram 
for different nanofilm thicknesses, where h and t are the 
normalized critical magnetic field h = iJc,||/77c,||(0) and 
the normalized temperature t = T/TffO), respectively. 
i7c,||(0) and rc(0) are the critical field at T = 0 and 
the critical temperature for iJ|| = 0. As presented in 
Fig. O for the nanofilm thickness d = 10 ML for which 
the orbital effect is negligibly small, the slope of h{t) 
at t = I is infinite and can be approximated by the for¬ 
mula dh/dtKi y^l/(I — Due to the orbital effect the 
slope h{t) at t = 1 becomes finite for thicker nanofilms 
(d = 20, 30 ML) and gradually decreases with increasing 
thickness. Similar behavior has been recently reported 
in the experiments with Pb nanofilm o^^di . The authors 
of Refs .^^di have argued that such dependence is caused 
by the boundary scattering, and therefore results from 
the roughness of the sample. As we have shown here the 
same behavior can be also induced by the orbital effect 
in the ultra clean film. 

The orbital effect in nanofilms influences also the or- 
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FIG. 5. (Color online) Spatially averaged energy gap A as 
a function of magnetic field and temperature T for the 
nanofilm thickness (a) d = 21 ML and (b) d = 23 ML. The 
value of the energy gap in each figure is normalized with re¬ 
spect to its maximum. 


der of the superconductor to normal metal phase transi¬ 
tion. It is well known that for the orbital limiting case, 
with no Pauli effect, the superconductor to normal metal 
transition induced by the magnetic field is of the second 
order. In contrast, in the Pauli limit the second order 
transition is suppressed with decreasing temperature and 
below t* = 0.56 it becomes of the first orde r^^d^ . The 
superconducting nanofllms are systems in which the or¬ 
bital and the paramagnetic effects are comparable while 
their relative importance can be controlled by changing 
the film thickness. By analyzing the free energy of the su¬ 
perconducting and the normal state we have determined 
the temperature t* at which the superconductor to nor¬ 
mal metal phase transition becomes of the first order for 
different nanofilm thickness. In Fig. [3 one can see that 
for ultrathin nanofilms, for which the paramagnetic effect 
is dominant, t* diverges to the value 0.56 predicted by 
Kazumi Maki in Refs 4^1^. If we increase the thickness, 
the temperature t* decreases which is directly related to 
the enhancement of the orbital effect. The extrapolation 
of t*{d) for t* = 0 gives d = 39 ML which is the limit 
above which the transition of the second order occurs in 



FIG. 6. (Golor online) Normalized critical magnetic field h = 
LIj, ii/TL^ ll(0) as a function of normalized temperature t = 
T/Tc{0) for different nanofilm thicknesses d. 



FIG. 7. (Color online) Temperature t* at which the supercon¬ 
ductor to normal metal transition becomes of the first order 
as a function of the nanofilm thickness d. 


the whole range of temperatures which corresponds to 
the orbital limiting case. 


IV. SUMMARY 


The superconductor to normal metal phase transi¬ 
tion driven by the in-plane magnetic field for Pb(III) 
nanofilms has been investigated in the framework of 
the BCS theory. We have shown that the orbital ef¬ 
fect suppresses the critical field as well as the ampli¬ 
tude of the critical field oscillations induced by the quan¬ 
tum confinement (shape-resonances). The analysis of the 
— Tc diagram allows to demonstrate that due to the 
orbital effect the slope iLc,|| — 2k at Tc(0) becomes finite 
and systematically decreases with increasing film thick- 
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ness. This result agrees with recent experiments for Pb 
nanofilmsHili. We have also analyzed the thermal effect 
and shown that the temperature t* at which the super¬ 
conductor to normal metal transition becomes of the first 
order reduces with increasing film thickness. For the ul- 
trathin nanofilms, in which the Pauli pair-breaking mech¬ 
anism is dominant, t* approaches to 0.56 in agreement 
with the theoretical prediction by Kazumi Mak i^^d^ . 

It is worth mentioning that our study take into ac¬ 
count only the electronic structure. However, in the 
nanoscale regime, the confinement affects not only the 
electronic spectrum but also the phononic degrees of free¬ 
dom which for nanofilms strongly deviates from that in 
the bulb^. The quantization of the phononic spectra 
in nanofilms and its influence on superconducting prop¬ 
erties have been considered in many papers^^— . More¬ 
over, the effect of the confinement on the electron-phonon 
coupling strength has been recently studied by Saniz et 
al. in Ref4^. In this paper— attractive electron-electron 
interaction has been derived with the use of the Green 
function approach beyond the contact potential approx¬ 
imation. It has been found that the increase of the criti¬ 
cal temperature observed in superconducting nanofilms is 
due to the increase of the number of phonon modes what 
results in the enhancement of the electron-phonon cou¬ 
pling. Since the modification of the phononic dispersion 
due to the confinement only slightly changes the super¬ 
conducting properties of nanofilms, the general results of 
the present paper remain valid. 

Nevertheless, there is another factor which can consid¬ 
erably affects the shape resonances in experiments. Note 
that in our study we assume the hard wall potential pro¬ 
file in the direction perpendicular to the film what refers 
to the case of the so-called free-standing nanofilms. How¬ 


ever, in experiments a thin film is never isolated but 
grows on a substrate (for Pb nanofilms usually on Si). 
Due to the interface effect the substrate layers strongly 
affect both the electronic structure and the phononic 
dispersion of the nanofilm. The role of the substrate 
on the shape resonances has been recently studied in 
Ref—. Taking into account a finite lifetime of the quan¬ 
tized states and modeling the substrate/thin-film inter¬ 
face by a more realistic finite step potential, the authors 
of Ref— showed that for the case of strong-coupling to 
the substrate the shape resonances are significantly sup¬ 
pressed in reference to the free-standing limit. However, 
for metallic superconductors (with long coherence length) 
the enhancement of the energy gap is so strong that the 
shape resonances should be experimentally observed de¬ 
spite the destructive influence of the substrate. 

Finally, we would like to point our that one should 
be careful when using the BCS theory in the description 
of the strongly-coupled superconductors such as Pb, as 
in some cases it can lead to overestimation of the size 
effect with respect to the behavior reported in the ex¬ 
periments^. The more appropriate description of the 
strongly-coupled superconductors requires the use of the 
Eliashberg theory. Nevertheless, the results from recent 
experiments^i^i^ d^d^ for Pb nanostructures are qualita¬ 
tively well described by the BCS model. 
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